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ABSTRACT 

The mass-loss rate of donor stars in cataclysmic variables (CVs) is of paramount 
importance in the evolution of short-period CVs. Observed donors are oversized in 
comparison with those of isolated single stars of the same mass, which is thought to be 
a consequence of the mass loss. Using the empirical mass-radius relation of CVs and 
the homologous approximation for changes in effective temperature T%, orbital period 
P, and luminosity of the donor with the stellar radius, we find the semi-empirical mass- 
loss rate M 2 of CVs as a function of P. The derived M 2 is at ~ 10" 9 - 5 -10" 10 M yr~ x 
and depends weakly on P when P > 90 min, while it declines very rapidly towards 
the minimum period when P < 90 min, emulating the P-T 2 relation. Due to strong 
deviation from thermal equilibrium caused by the mass loss, the semi-empirical M 2 
is significantly different from, and has a less-pronounced turnaround behavior with P 
than suggested by previous numerical models. The semi-empirical P-M 2 relation is 
consistent with the angular momentum loss due to gravitational wave emission, and 
strongly suggests that CV secondaries with 0.075 M Q < M 2 < 0.2 M are less than 2 
Gyrs old. When applied to selected eclipsing CVs, our semi-empirical mass-loss rates 
are in good agreement with the accretion rates derived from the effective temperatures 
T\ of white dwarfs, suggesting that M 2 can be used to reliably infer T 2 from T\. Based 
on the semi-empirical M 2 , SDSS 1501 and 1433 systems that were previously identified 
as post-bounce CVs have yet to reach the minimal period. 

Subject headings: binaries: close — binaries: eclipsing — stars: novae — cataclysmic 
variables — stars: dwarf novae — stars: fundamental parameters — stars: low-mass — 
stars: mass loss — white dwarfs 



1. Introduction 



Cataclysmic variables (CVs) are semi-detached binaries with a white dwarf (WD) primary and 
a low-mass companion which transfers mass to a WD through Roche-lobe overflow. The observed 
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orbital period distribution of CVs has three important features: a sharp short period cut-off at the 
orbital period of P m i n — 76.2 min, referred to as the "minimum peri od", the perio d gap, a dearth of 
systems with periods between -P gE ^> — 129 min and -P ga Jp — 192 min (jKnigg e 2006), and the "perio d 
minimum spike", a significant accumulation of systems at P ~ 80-86 min (jGansicke et al.ll2009l ). 
For CVs with period above -P gE ^, the dominant angular momentum loss mechanism is thought 
to be the magnetic braking, while the gravitational radiation may be responsible for the angular 
momentum loss for CVs below P~. In this work, we focus on CVs only with P < P~. 

gap * J gap 

A characteristic feature of CVs is the mass transfer from a secondary to a WD. Since the 
gas stream carries substantial angular momentum, it forms an accretion disk around the WD. 
The reaction of a CV on the mass transfer can be qualitatively understood with the help of the 
dimensionless parameter 

T = — , (i) 

where tm and tk h denote the mass-loss and Kelvin-Helmholtz timescales of the donor, respectively 



(jPaczynskil 1 1 98 ll ) . At the initial phase of mass transfer, t> 1, so that the donor remains close to 



thermal equilibrium and thus follows the mass-radius relation (MRR) of main-sequence stars 

R 2 oc M| E , 



(2) 



with the effective mass-radius exponent £e close to 0.64^q'q2 ( Kniggd 120061 ). Here, the subscript 
"2" denotes the quantities of CV donors. Both R 2 and P of the CV decrease as the donor keeps 
losing mass during this early phase of CV evolution. 

As the mass transfer continues, the binary shrinks due to the angular momentum loss, reducing 
r. When tm and tkh become comparable to each other (r ~ 1), the donor becomes out of thermal 
equilibrium. Since the surface luminosity is no long er in balanc e with the nuclear luminosity, the 
donor expands, causing £e to decrease with ti me ([Kind Il988l ), Using the Hayashi theory and 
homologous approximation, IStehle et al.l (jl996l ) studied the reaction of the donor to mass loss 
and found that R 2 increases with decreasing r. Observations confirm this theoretical prediction, 
showing that the ra dii of donors in short-period CVs a re ~ 5-50% larger than those of isolated 



main-sequence stars (jPatterson et al 



2005 



Kniggd 12006J). 



The effective mass-radius exponent £e becomes smaller as M 2 decreases towards the brawn 
dwarfs masses. The change in the MRR of donors slows down the temporal change in the orbital 
period as the donor continues los ing mass. The orbital period reaches a minimum value when 
£ E k, 1/3 at which dP/ dM 2 = (jRappaport et al.lll982l ). As r decreases further, £e declines to 



0.21^q'5q (|Kniggdl2006l ). This in turn leads to an increase in P. Such a turnaround of the orbital 



period at Ptum is referred to as the "period bounce" . 

The evolutionary scenario described above suggests that the mass-loss rates M 2 (and the 
related thermal relaxation) of donors control the evolution of CVs. It also implies that the observed 
distribution of CVs versus the orbital period may be a consequence of a certain relation between 
M 2 and P. Since M 2 (P) not only provides a way to estimate the ages of CVs but also can be used 
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to check the dominant mechanism for angular momentum loss, finding a reliable P-M2 relationship 
is of great importance in understanding CV evolution. 

There have been numerous studies to deri ve the P—M?, relation; these can be categorized 



Rappaport et al.l Il982l: 



into two groups: purely theoretical w ork (e.g , 
Kolb fc Bara"ffelli999l: Ho well et al I hooih and semi-empirical work (e.g 



Ritter 



Patterson 



198£ 



i Jem 



1988 



19841 : llTrban fc Sion 



20061 : Townslev k, Gansickell2009l ). The first, theoreti cal studies ma de use of stellar evolution codes 
combined with the analytical mass-transfer models of lRitterl (j 19881 ) or others. The common results 
of these numerical experiments are that (1) the orbital peri od exhibits the turnaround behavior at 



turn 



that is about ~ 10 % smaller than the observed P m \ n (jRappaport et al.lll982l ; iKolb Sz Baraffe 



1999l : lHowell et alJboOll ) and (2) compared to isolated stars, the secondaries in short-period CVs are 
oversized due to the mass loss but only by ~ 5-20%, which is about a factor ~ 1.5-2 smaller than 
the observed CV sizes. Since the bloating of the donor depends on the efficiency of the mass loss, 
these results suggest either the numerical models unde restimate above P V m, or the efficiency o f 



the angular momentum loss is higher than expected (jKolb fc Baraffdll999l : 



Renvoize et al 



Desp i te the attempts to include the effects of t idal and rotational perturbations (jRezzolla et. al 



2002). 



2001 



Renvoize et al 



2002: 



Kolb &; Baraffdll999l ). additional angular momentum loss via magne tic 



braking (lAndronov et al.ll2003l ). accretion disk winds and magnetic propeller (see lBarker 1120031 for 
review), etc., clear explanations for the discrepancies between Ptum and -P m i n and between the 
empirical and numerical MRRs have yet to be found. 

The second, semi-empirical approach uses, for instance , the effective temperatures of WDs 



dTownslev fc Bildstenlhool . [2004J; iTownsley &: Gansickdl2009i ) or accretion disk luminosities ([Patterson! 



19841 ) to infer the accretion rates M\. The physical basis for the for mer is that WDs in CVs 
are observed to be hotter than expected from their ages (jSionl Il995l ). and such "overheating" 
of WDs may be a result o f the compressional heating due to the gas piled up on their surfaces 
(ITownsley &; Bildstenl 120031 ) . For the latter method, the accretion luminosities can give direct in- 
formation on Mi if the disk properties are well constrained. When applied to dwarf novae (DNe) 
that accrete sporadically and have thermally unstable disks, M\ calculated from the two methods 
have a different scatter for given P. The mass-transfer rates based on the accretion luminosities 
probably measure an accretion rate averaged over short ti me intervals, e.g., a few deca des, which 
makes them quite sensitive to short-time disk variability (jTownsley Gansickd l2009i ) . This in- 
evitably leads to a high dispersion in Mi estimated from the disk luminosities. On the other hand, 
the effective temperatures of WDs are likely to give a relatively long-term averaged M\ , although 
this method may suffer from uncertainties related to unknown properties such as masse s and radii 
of W Ds as well as some auxiliary assumptions on their structures and boundary layers ([Patterson 
2009T ). 



In this paper we construct another semi-empirical relation between the mass-loss rate and the 
orbital period, and compare the result with those of the purely t heoretical and the semi-empirical 
approaches mentioned above. We start from the empirical MRR of lKniggd (120061 ) for superhumping 
CVs, and apply the homologous approximation in order to calculate the dependencies of the effective 
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temperature and gravothermal luminosity upon the bloating of the donor. Assuming that the 
observed MRR of CV donors is a consequence of the mass loss and thermal relaxation processes, 
we calculate the mass-loss rates. We will show that the resulting mass-loss rates are broadly in 
line with those on the basis of t he effective tem p eratu res of WDs, while different from those of 



the purely theoretical models of iKolb fe Baraffd (|1999l ). We will also show that our results are 



consistent with the angular momentum loss of CVs due primarily to the emissi on of gravitat ional 



waves. Since radii of CV donors likely trace the long-term evolution of M2 (IKniggd 120061 ) and 
depend only weakly on details of the mass-transfer processes, the resulting mass-loss rates from the 
empirical MRR are less model-dependent and less sensitive to the short-term effects than the other 
semi-empirical estimates that rely on accretion activities. 

The organization of this paper is as follows. In £|2j we introduce the theoretical MRR for 
isolated single stars and the empirical MRR for CVs donors, and define the bloating factor. We 
then describe a method to find the effective temperature and gravothermal luminosity of donors 
by using the bloating factor. In Sj3l we calculate the semi-empirical mass-loss rate, construct the 
mass-loss history to estimate the CV ages, and show that the derived mass-loss rate is consistent 
with CV evolution by angular momentum loss via the emission of gravitational waves. In £jH we 
apply our results to selected eclipsing short-period CVs and show that the semi-empirical mass-loss 
rate is consistent with the accretion rates inferred from the effective temperatures of WDs. We also 
suggest a way to estimate the effective temperatures of donors assuming that mass is conserved. 
Finally, we summarize our results in Sj5j 



2. Mass-radius and Period-effective Temperature Relations 



We consider a CV consisting of a WD (or accretor) with mass M\ and a secondary (or donor) 
with mass M2, radius R2, and effective temperature T?, . We assume that the donor fills its critica l 
equipotential surface. In accordance with the results of lPatterson et al.l (120051 ) and IKniggd ([2006), 
we assume M\ = 0.75 M throughout this paper. We use an additional subscript "0" to denote the 
quantities of isolated, single stars in hydrostatic and thermal equilibrium to which CV secondaries 
will be compared. In this section, we use the empirical relationship between masses and radii of CV 
secondaries to calculate the bloating factor in comparison with the single stars of the same mass. 
We then treat the bloating factor as perturbations to the properties of single stars in order to find 
the orbital period, effective temperature, and gravothermal luminosity of CVs. 



2.1. Mass-radius Relations 



As the MRRs of single stars that we use as the refe rences t o thos e of CV secondaries, we take 
composite relations constructed in a manner similar to IKniggd (|2006l ): the BCAH9 8 isochrone of 
Baraffe et all (|l998h for stars with 0.075 M Q < M 2 < 0.3 M Q , the DUSTY isochrone of lBaraffe et al 
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(12002) for 6 < M 2 < 0.075 M (1500 < T 2 < 2000 K), and the COND isochrone of 
Baraffe et all (|2003l ) for M 2 < 0.06 M . 



There are quite large uncertainties in the estimation of CV ages. Period bouncers are thought of 
as a final stage o f CV evolution, suggesting that they are olde r than their "parent" CVs (see, e.g., 
Paczynskil Il98ll : iRappaport et al.l Il982l : iKolb &: Baraffd 1 1999I ). Although numerical experiments 
predict that CVs cross the perio d gap at ~ 1 Gyr and reach Ptum m ~ 4 Gyr after the onset of 
the mass loss below the gap (e.g., IKolb k Baraffd Il999l : iHoweU et allboOlh . temporal evolution of 
CVs depends strongly on the mass-loss rate that is not well constrained observationally. To cover 
a wide range of CV ages, we adopt two composite isochrones with 1-Gyr and 5-Gyr ages (hereafter 
Bl and B5 sequences, respectively )0. Figure Q] plots the composite i?2,o of single stars as functions 
of M 2 for the Bl and B5 sequences, as dashed and solid lines, respectively. Both sequences give 
almost the same radii above 0.075 Mq, while below this point the stellar radii in the Bl sequence 
are about ~ 10% larger than in the B5 sequence. 



As the MRR of CV secondaries, we take the relation derived by iKniggd (|2006l ) from the 
observed properties of superhumping CVs: 



R2 

Rtn 



0.1101 



M 2 
0.063 M Q 



0.21 



0.063 M < M 2 , 

0.063 M < M 2 < 0.2 M , 

0.2 M < M 2 . 



(3) 



This empirical fit assumes Mi = 0.75 M and incorporates itself the additional empirical constraints 
that reproduce the observed locations of the period gap and the period minimum, such that -P gap 
occurs at M 2 = 0.2 ± 0.02 M and P min at M 2 = 0.063 ± 0.009 The intrinsic errors in R 2 
for this empirical MRR are estimated to be about 2-3% (jKniggd 120061 ) ; we take allowance for 4% 



uncertainties in R 2 throughout this work. We note a caveat that our choice of the errors may 
underestimate real values below 0.075 M where only a few systems are observed and thus there 
is practically no calibrator. In addition, the masses of WDs may differ from the assumed mass of 
0.75 M , which also can change R 2 . Equation ([3]) is plotted in Figure [T] as dotted line with the 
shaded region denoting the associated uncertainties in R 2 . Columns (1) and (2) in Table[T]give the 
M 2 -R 2 relation. 

For comparison, Figure [TJ also plots as a dot-dashed line the MRR of CVs with 0.03 M < 
M 2 < 0.21 M from the numerical models of lKolb &: Baraffd (|1999l ). which clearly has smaller radii 
than the empirical relation, typically by ~ 5-10%. This indicates that the deviation of real donors 
from thermal equilibrium is larger than that expected from the numerical models. 



1 The composite isochrones we adopt are slightly different from those in Kniggel ( 2006) who took the 5-Gyr is ochrone 
for do nors with M > 0.075 Mq and 1-Gyr isoch rone for M < 0.075 Mr . The choice of the 1-Gy isochrone in Kniggel 
(2006) was to represent the numerical results of Kolb fc Baraffe 1 1999 



- 6 - 



0.40 




0.07 I , , I i i , 

0.03 0.1 0.3 

M 2 /M 



Fig. 1.— Ma ss-radius relations of CV donor stars. The dotted line plots the empirical MRR from 
Kniggd (|2006l ) with the 4% uncertainties indicated by the shaded region. The dashed and solid lines 



represent the unperturbed MRR of singles stars for the Bl and B5 sequences, re spectively. The 
dot-da shed line is for the MRR of evolving donors from the numerical models of iKolb &: Baraffe 
(|1999l ). The thin vertical lines mark the boundaries of the adopted isochrones: BCAH98 for 
0.075 < M 2 < 0.3 M ; DUSTY for 0.06 < M 2 < 0.075 M ; COND for M 2 < 0.06 M . 
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2.2. Bloating Factor 



We define the bloating factor 



a 



R2.0 ' 



(4) 



as the relative size of a CV secondary to a single star with the same mass. Figure [2] plots a as 
dashed and solid lines for the Bl and B5 sequences, respectively, with the shaded and hatched 
regions representing the corresponding uncertainties. The discontinuity of a at M 2 = 0.2 M Q 
corresponds to the period gap. Since P oc a 3//2 (see below), the ratio of the bloating factors above 
( a gap) and below (ag~ ap ) the gap satisfies «+, /a~ ap = (P g & p / Pg&p) 2 ^ = 1-3, in good agreement with 
the width of the gap. At the upper edge of the gap, the donors detach from their critical lobes and 
the mass loss stops. Since the donors are thought to be oversized due to the mass loss, they shrink 
within the gap. When the donors come out the lower edge of the gap, they resume the mass loss 
and become oversized again. 

Figure [2] shows that the bloating factor increases with the decreasing donor mass for M 2 > 
0.075 M & and then decreases with decreasing M 2 below this point. This can be explained qual- 
itatively as follows. Since stars with M 2 < 0.3 M © are fully connective an d isentropic, they are 
well described by polytropes with index n = 3/2 (jBurrows fc Leibertl Il993l ) . The MRR for these 
polytropes 

R 2 oc KM~ 1/3 , (5) 



with K being the pressure constant (|Chandrasekharlll939l ). then gives 



a 



K 



(6) 



The mass loss is capable of driving the donor out of thermal equilibrium. This raises the gravother- 
mal luminosity, L g = L 2 — L nuc , where L nuc is the nuclear luminosity of the donor and L 2 is the 
surface luminosity. For fully convective stars, the change of K depends on L g and tkh as 



K 
K 



1 L e 



(X 



(|King 



1988 



tkh L 2 

Ritteiill996l ). Combining equations ((5J), ©, and ((7J), one obtains 



(7) 



a 1 L s 1 /. 1 

- oc -r + — Ueo+ ~ 

Oi T K H ■L'2 Tm \ 3 



(8) 



where £eo is the mass-radius exponent of unperturbed, singe stars. At the lower edge of the period 
gap, a donor remains close to thermal equilibrium so that L g ~ and £eo ~ !• I n this case, 
the second term in the right-hand side of equation (JSj) dominates and a/ a > 0, indicating that a 
increases as the mass loss continues. When the donor reaches the minimum hydrogen burning mass 
(MHBM) ~ 0.075 M©, the nuclear burning stops and donors are born as brown dwarfs, so that 
Lg/L 2 = 1. In the brown dwarf regime, on the other hand, the internal pressure is dominated by 
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Fig. 2. — Bloating factors a of CV donors as functions of their mass M%. The dashed and solid lines 
are for the Bl and B5 sequences, respectively, with the shaded and hatched regions corresponding 
to the uncertainties. The dotted line shows the bloating factor of Roche-lobe filling polytropes 
with n = 3/2 due to the isentropic expansion caused by tidal and rotational perturbations, when 
its companion has a mass of M\ =0.75 M . 



the d egeneracy pressure with K depending mainly on the chemical composition ([Burrows Leibert 
19931 ). In this case, £eo tends to —1/3 and the first term in the right-hand side of equation (J5J 
dominates, resulting in a/ a < 0. For brown-dwarf donors, the change in the equation of state 
partially compensates for the increasing tendency of the bloating factor due to the mass loss. 

In addition to the mass loss, tidal and rotational perturbations can also change the donor size, 
depending on the polytropic index and the mass ratio. We assume that the applied perturbations do 
not change the specific entropy, which is reasonable since the stellar distortions occur in dynamical 
timescales. In Appendix [A] we calculate the enlargement factor ai sen of a Roche-lobe filling donor 
as a function of n for q = 0.1. The resulting ai sen for n = 3/2 polytropes against M<i is plotted in 
Figure[2]as a dotted line. Note that the expansion of low-mass donors due to the ti dal and rotational 
distortions is ~ 4-5%, insensitive to the mass ratio for q ~ 0.04-0.4 (see also, e.g.. lUryu fe Eriguchi 



1999 



Renvoize et al.1 120021 : ISirotkin fc Kimll2009l ). It can be seen that while ai se n is, in general, 
considerably smaller than the observed bloating factor for stars with 0.04 Mq < M<i < 0.15 M , it 
is comparable to a when M% ~ 0.2 M corresponding to P g ^ p or when M2 < 0.04 Mq as donors 
approach P m i n . This indicates that the effect of the isentropic expansion of donors on the bloating 
factor is non-negligible for CVs that just resume mass loss or are below the minimum period. 



2.3. Orbital Period and Effective Temperature 



For a given MRR, the orbital period can be determined by the Roche lobe geometry and 
Kepler's law as 



P 2 



3vr 



Gpp(l + qY (9) 

where p = 3M2/ '^irR^ is the mean density of the donor, / = (2/3 4 / 3 )g 1//3 /(l + g) 1 / 3 i s the relativ e 
size of the donor to the orbital separation, and q = M2/M1 is the mass ratio (e.g., lKopallll972l ). 
Equation implies that the orbital period depends on the bloating of the donor as 



P 



Pn« 3/2 . 



(10) 



Clearly, the orbital period increases with the radius of a secondary. 



Since donors below the period gap are fully convective, their lumino sities are determined by 
the conditions at the very outermost layers (jKippenhahn &: Weigertlll990l ). so that one can neglect 
the influence of the distortions on the central parts of a donor. To obtain the properties of the 
photospheric layers, we use the standard Hayashi theory. For the opacity law of the form 



k = K P^T°, 



(11) 



with the gas pressure P g , the ratio of the effective temperatures satisfies 
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where 7 = (3a — l)/(5a + 26 + 5) (e.g, iKippenhahn k, Weigertlll99Cl . see also IStehle et al.lll996l ). 
For the realistic opacity laws applicable for low-mass main-sequence stars, 7 has a very small 
value. If the opacity is governed by H~, for example, a « 0.5 and b « 4.5, leading to 7 w 0.03. 
If the molecular absorption instead is a dominant opacity source (c w 1 and b « 0), 7 « 0.2. 
Therefore, donors in short-period CVs have approximately the same effective temperatures as the 
corresponding isolated single stars of identical masses (e.g.. IStehle et al.l ll996). 



To construct the semi-empirical P—T2 relations, we first use the theoretical M2—T2 relation 
for isolated stars and calculate Po from equation ([9]) by fixing Mi = 0.75 M . We then use the 
empirical MRR and the bloating factor of donors to calculate P and T 2 from equations (|10p and 
(|12p . The Pq— T 2) o relations are plotted in Figure [3] as thin lines, while the thick lines are for the 
semi-empirical P-T 2 relations. The dashed and solid lines correspond to the Bl and B5 sequences, 
respectively. Again, the shaded and hatched regions give the corresponding uncertainties. Columns 
(3), (6), and (11) in Table Q] list the values of P and T 2 against M2 for both sequences. For a given 
mass, the expansion of a donor increases the orbital period while the effective temperature keeps 
almost unchanged. The largest difference between P and Po occurs at T2 ~ 2000 — 2200 K, 
corresponding to M2 ~ 0.07 5 where th e bloating factor is maximal. The dot-dashed line plots 



the semi-empirical resul t of iKniggd (120061 ). Because of the difference in the adopted ages of the 
isochrones, the result of kniggel feood ) for M 2 > 0.075 M (with T 2 > 2000 K) and M 2 < 0.06 M Q 
(with T2 < 1800 K) is in agreement with our results for the B5 and Bl sequences, respectively. 
The donors at P min have mass M 2 ~ 0.054 - 0.072 M Q , corresponding to T 2 ~ 1000-1800 K and 
~ 1600-2100 K for the B5 and Bl sequences, respectively. The right choice of the stellar age is 
quite import ant in evaluating the m ass-loss rate near Pjxxin since it depends rather sensitively on T2 
as M 2 cx T 2 4 (|Rappaport et al.lll982l ). 



2.4. Luminosities and Timescales 

Combining equation (|12p with the Stefan-Boltzmann law 

L 2 = 4vra^T 2 4 , (13) 

where a is Stefan-Boltzmann constant, one can show that the change in the stellar luminosity is 
given by 

22q+4b+6 

-^2/-^2,0 = « 5 a + 2i >+5 . (14) 

Although the effective temperature does not change much with the bloating factor, the additional 
o? term in the Stefan-Boltzmann law increases L 2 /L 2; o significantly. Columns (5) and (10) in Table 
[Ogive the values of L 2 for the Bl and B5 sequences, respectively. 

The nuclear luminosity L nuc is determined by the energy produced in the central parts of a 
donor, so that L nnc ~ M 2 e c , where e c is the thermonuclear reaction rate given by 

e c oc T^- 1 , (15) 
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Fig. 3. — Dependence of the effective temperature T2 of a donor on the orbital period P. The 
thick dashed and solid lines plot the semi-empirical P-T2 relations for the Bl and B5 sequences, 
respectively, with the dashed and hatched regions denoting the corresponding uncertainties. The 
thin dashed and solid lines give the unperturbed Pn~^2,p relat ions for the Bl and B5 sequences, 
respectively. The semi-empirical P-T2 relation of iKniggd (|2006l ) is plotted as the dot-dashed line. 
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with T c and p c denoting the central temperature and density, respectively. Since T c oc oc a 

for a monatomic ideal gas, the change in the nuclear luminosity due to the bloating factor is 

Lvmt a 3 ". (16) 



L 



nuc,0 



Theo retical models for stellar structure and evolution give v ~ 6 and rj ~ 2 (jBurrows fc Leibert 
19931 ) for T c ~ 3 x 10 6 K and p c ~ 150 g cm -3 , typical for the centers of low-mass stars. 



Since £2,0 = L nuC £ for isolated stars in thermal equilibrium, the gravothermal luminosity 
satisfies 

rl = 1 _ nuc (17) 

L2 £nuc,0 ^2 

Thermonuclear burning stops in CV secondaries once they reach the MHBM ~ 0.075 M . At 
this point, donors are effectively born as brown dwarfs and the gravitational collapse becomes the 
dominant source of the stellar radiation. That is, L g /L2 = 1 for M2 < 0.075 M©. 

Figure H] plots L^jL^ versus M2 for the Bl and B5 sequences as dashed and solid lines, respec- 
tively, with the uncertainties denoted by the shaded and hatched regions. Figures [2] and H] together 
imply that at any given mass the contribution of the gravothermal luminosity increases with the 
bloating factor. If R2 increases by 20%, more than 80% of the donor luminosity is generated by 
the release of thermal and gravitational energies. 



3. Mass-loss Rate 



Finding the theoretical mass-loss rates of CV donors is a daunting task since it involves many 
complicated processes including thermal and adiabatic relaxations, nuclear evolution, effects of tidal 
and rotational distortions, orbital angular momentum loss, etc., al l of wh ich may affect the internal 
structure significantly. Although the stationary model of iRitterl (|1988l ) has been used in various 
theoretical models of CV evolution, it relies on the parameters describing the donor atmosphere 
in the vicinity of the inner Lagrangian point, which are poorly known. By relating the mass-loss 
rate with angular momentum loss, on the other hand, the Ritter model allows to check whether a 
proposed angular momentum loss mechanism is indeed responsible for CV evolution. As mentioned 
in ffl t he nu merical calculations that incorporate the stellar evolution and the mass-loss model of 
Ritterl (119881 ) produced MRRs that differ from the empirical MRR. 



3.1. Semi-empirical Mass-loss Rate 



Instead of employing the Ritter model to calculate the mass-loss rate, we seek the M2 that is 
required to reproduce the empirical MRR. We make use of the well-known relation between the 
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Fig. 4. — Semi-empirical ratio of the gravothermal luminosity L g to the surface luminosity L2 as 
a function of the donor mass. The dashed and solid lines correspond to the Bl and B5 sequences, 
respectively, with the shaded and hatched regions denoting the uncertainties. For donors below the 
MHBM, all radiated energy comes from gravothermal collapse (L g /L2 = 1). 
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effective mass-radius exponent £e of the donor and the mass-loss timescale: 

t,E = t,S + T 77, 
o Li 



(18) 



(|Ritterlll996l ). where £s ~ —1/3 is the adiabatic mass-radius exponent for fully convective stars. 
Equation (I18p dictates how the MRR should behave in response to the mass-loss and thermal 
relaxation processes. When the mass transfer just commences, L g ~ and the donor response is 
adiabatic, leading to ~ £s- Since the ensuing deviation from the adiabatic evol ution depend s 
on how the donor reacts thermally on the mass loss, £e increases linearly with r (jRitterl 11996), 
Combining equations ([9]), (I13p and f 1 1 8 [) with the definitions of the Kelvin- Helmholtz timescale 
(tkh = GM% /(R2L2)) and mass-loss timescale (tm = M%j \ M2 I), one obtains 

(Ce + 1/3) 243 7T L 2 2 V ; 

1/3 corre sponding to the period bounce, equation (fl~9|) recovers equation (34) of 



When £p, 



Rappaport et al.l (jl982l ) for the mass-loss rate at the minimum period. 



In Figure[5]we plot the semi-empirical P-M2 relations calculated from (|19p for the 1-Gyr and 5- 
Gyr isochrones as dashed and solid lines in red color, respectively. The corresponding uncertainties 
are given by the shaded and hatched regions. Columns (7) and (12) of Table CD list the values 
of Mi- The semi-empirical mass- loss rate is essentially non-stationary. For 129 min > P > 90 
min, M2 is in the range of ~ 10~ 9,5 - 10 -10 M yr _1 and decreases relatively slowly as M2 oc P 3 . 
Below 90 minutes, M2 sharply declines towards P m \ n as a result of a strong dependence of M2 
upon T2, emulating the P-T2 relation shown in Figure [3l Since the accretion rate decreases with 
M2, this rapid decline of M2 is consistent with th e observed low accreti on activities of CVs within 
the period spike identified in the SDSS samples (jGansicke et al.ll2009l ). Note that the mass- loss 
rates for the Bl and B5 sequence are almost similar to each other for masses above the MHBM, 
suggesting that these are perhaps permanent values. In the brown dwarf regime, on the other hand, 
the ambiguity in the adopted stellar age leads to a large difference in the mass-loss rate, suggesting 
that it is necessary to use the appropriate isochrones for an accurate estimation of M2 for donors 
with M 2 < 0.075 M . 

For c omparison, Figur e [5] als o plots as a dot-dashed line the P-M2 relation from the numerical 
models of iKolb &: Baraffd (|1999l ). which shows a pronounced "boomerang" shape such that M2 
is nearly stationary above Ptum with M2 oc P ' 5 . It also has M2 significantly smaller for M2 > 
0.075 Mq than our semi-empirical results. To see what causes the difference between our semi- 



empirical P-M2 relation and the n umerical result of IKolb &: Baraffd (|1999l ). we calculate M2 from 
equation ([19]) by using the MRR of IKolb & Baraffd* |l999h together with the 5-Gyr isochrone. The 
resulting P-M2 relation is pl otted in Figure [5 ] as a dotted line, which is overall in good agreement 
with the numerical results of lKolb &: Baraffd (|1999l ). This suggests that smaller radii of donors in 
their work are a prim ary reason for the discrepancies between the results of the current work and 
Kolb fc Baraffd (|1999l ). This also proves to some extent that our approach provides a reasonable 
way to estimate the mass-loss rate. 



60 70 80 90 100 110 120 130 

P (min) 



Fig. 5. — Dependence of the mass-loss rate M2 of a donor on the orbital period P. The dashed 
and solid lines in red give M2 for the Bl and B5 sequences, respectively, with the dashed and 
ha tched regions denoting the uncertainties. The dot-dashed line represents the numerical result 
of iKolb Baraffd ()1999l ). while the dotted line plots the mass- loss rate obtained from equation 
(|19p with the numerical MRR of lKolb Baraffd (|l999l ) for evolving donors. The dashed and solid 
lines in green den ote Mob., for the Bl and B5 sequences, respectively, from the stationary mass- 
transfer model of iRitterl (|1988l ) when the angular momentum loss is dominated by gravitational 



wave emission. 
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From the semi-empirical M2, we calculate r for the Bl and B5 sequences, which are listed 
in Columns (8) and (13) of Table [fl respectively. Note that donors in CVs below P g ~ p are out of 
thermal equilibrium with r < 1, which is the main cause of the donor expansion. Equations (|17p 
and (|19p imply that for given Mi and £e 5 the mass-loss rate increases more than c? due mainly 
to changes in the orbital period and gravothermal luminosity. The resulting semi- empirical P-M? 
relatio n still has a boomerang shape, which is much less dramatic than that of iKolb &z Baraffe 
(|1999l ) due to the donor bloating and the associated thermal relaxation. Nevertheless, it is very 
difficult to detect the boomerang shape observationally because of the sharp decline of T2 and M2 
below 90 min. 



Once the P-M2 and P-M2 relations are found, it is straightforward to calculate the temporal 
change of the donor mass by solving dM2/dt = We assume that donors are 1 Gyr old and 
have a fixed initial mass of .2 Mm when they e merge from the period gap. Since CVs with P < 2 
h are predominantly DNe (jRitter &: Koibll2003l ). we further assume that the masses of WDs are 
constant at M\ = 0.75 M throughout CV evolution. Figure [6] plots as dashed and solid lines 
the evolution of donor mass and orbital period calculated from the semi-empirical M2 for the Bl 
and B5 sequences, respectively, as functions of time t elapsed si nce P~ (i.e., tp = 1 Gyr). The 



dot-dashed lines show M2{t) and P(t) from the mass-loss rate of lKolb Baraffd (|1999l ) for which 



Mi = 0.6 Mq and M2 = 0.21 M m initially at P^ ari - Th e changes in the donor mass and orbital 



period in the numerical model of IKolb &: Baraffd (|1999l ) are in general slower because of a lower 
mass-loss rate. Our semi-empirical results indicate that M2 decreases from 0.2 Mq to 0.075 Mq in 
~ 1 Gyr and that the minimum period is reached in ~ 1-2 Gyr after -P ga T), with out no essential 



depen dence on the adopted stellar age. On the other hand, the numerical model of lKolb Sz Baraffe 
(|1999l ) predicts that the same decrement of M2 occurs in ~ 3.5 Gyr, while it takes the donors ~ 4 
Gyr to reach P m i n . Although uncertainties in the stellar ages make it difficult to follow M2{t) and 
P(t) accurately below the MHBM, the close agreement between the results based on the Bl and B5 
sequences suggests that the ages of CV secondaries with 0.075 M Q < M2 < 0.2 M Q are less than 2 
Gyr. 



3.2. Consistency check 



Mas s transfer in CVs below the period gap is considered to be driven by grav itational wave ra- 



diation ( Paczynski 



1981 



Patterson 



1984; 



Kolb fc Baraffe 



1999 



Howell et al.ll200ll ). While equation 



(|19p is useful to obtain M2 from the empirical MRR independent of specifics of mass transfer, it 
does not provide any clue regarding what dri yes mas s tran sfer in CVs. In this subsection, we make 
use of the stationary mass-transfer model of iRitterl (119881 ) and check whether the semi-empirical 
M2 calculated in E j3.ll is consistent with the angular momentum loss via emission of gravitational 
waves. 



Ritteri ()1988l ) presented a simple analytical model in which the mass transfer from a donor is 
treated as a stationary isothermal, subsonic flow of gas through the inner Lagrangian point. Taking 
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Fig. 6. — Temporal evolution of (a) the donor mass Mi and (b) the orbital period P due to mass 
loss. Note that t denotes the time measured from the lower bound -P ga ^, of the period gap where 
the donor has an age of to = 1 Gyr. The dashed and solid lines are for the Bl and B5 sequences, 
respectively, with the shaded and hatched regions denoting the u ncertaities. The dot-dashed lines 
show the results from the mass- loss rate of lKolb fc Baraffd (119991 ) . 
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allowance for the finite scale height of the donor photosp here and cons idering the temporal changes 
in the volume radii of the donor and the Roche lobe, iRitterl (|1988l ) showed that the stationary 
mass-transfer rate is given by 



M- 



2,R - 77 Vt oT^ > ( 20 ) 

where rj is the angular momentum loss timescale, £r is the mass-radius exponent of the critical 
Roche lobe. Note that we neglect the nuclear-timescale term in equation (|20p since it is much 
longer than the other timescales involved. In the case of DNe for which the mass transfer is 
non-conservative, £r, is given by 



7U 1 



+ 



+ 2q 



(21) 



:! Jl + q 

(jSoberman et al.l 119971 ). If the gravitational wave emission is the dominant mechanism for the 
shrinkage of the orbital radius, tj can be set equal to tq defined as 



/ P\ 8/s (M 1 + M 2 ) 1 / 3 



32 G 5 /3 \2tt) 



M X M 2 



(22) 



(jLandau & Lifshit j|l975l ). 



We use the empirical MRR of iKniggd (|2006l ) and solve equation (|20p with rj = tq for M 2 ^. 
The resulting P-M 2 ^ relations are plotted in Figure [5] as dashed and solid lines in green color for 
the Bl and B5 sequences, respectively. The adopted stellar ages make a small difference in M23, 
since the first term in the square brackets of equation f)20f) is insensitive to the isochrones when 
M 2 > MHBM, while it becomes much smaller than the second term when M 2 < MHBM. Figure [5] 
shows that M 2 r agrees, within errorbars, with M 2 for both Bl and B5 sequences before reaching 
the minimum period. This verifies that our semi-empirical mass-loss rate is consistent with the 
angular momentum loss due to emission of gravitational waves, at least for pre-bouncers. Beyond 
P m in, M 2 for the Bl sequence appears to match M23 better than the B5 sequence, suggesting that 
the CV ages are presumably less than 5 Gyr. 



4. Discussion 



Since the mass-loss rate cannot be observed directly, we are only able to compare our semi- 
empirical results with the accretion rates estim ated by means o f other i ndicators. For this pu rpose, 
we select 10 short-per i od, ec lipsing C Vs from iPatterson et al.l ()2005l ) , iLittlefair et al.l (J2008J) , and 
Townsley &: Gansickd (J2009J). Table [2] represents a compilation of our sample CVs, with the pa- 
rameters in Columns (2) to (7) taken from the referred papers. To our knowledge, these are the 
only eclipsers with P below 130 min, for which the reliable data for the binary parameters such as 
P, Mi, M 2 , and T\ are all available. To infer the mass-accretion rates for the sample CVs, we use 
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where T\ and M\ are the effective temperature a nd m ass of a WD and (Mi) is the time-averaged 
accretion rate, taken from iTownslev Gansickd ([20090 . Since the gravitational energy released by 
the accreted material onto the surface of a WD is deposited at shallow depths and soon radiated 
away, the quiescent surface luminosity should be proportional to the mass-accretion rate, leading 



1/4 



to T\ oc (Mi) (e.g.. ITownslev Sz Bildsten 



2003 



2004 ). Filled circles with errorbars in Figure 0a 



plot (M) for the sample CVs calculated from equation ([23]) ; these are also listed in Column (8) 
of Table [2j It can be seen that (Mi) is broadly in line with our semi-empirical M2 for both Bl 
and B5 sequence. This suggests that our semi-empirical P- M 2 relation is realistic, p rovid ed that 
—Mi = (Mi) in these CVs during accretion disk quiescence. ITownslev fc Gansickd (|2009i ) noted 
that (Mi) significantly declines towards -P m i n for P < 84 min. Since the accretion rate decreases 
with decreasing M2, this downward tendency of (Mi) with P is in agreement with the behavior of 
the semi-empirical M2. 

Donors in short-period CVs are very faint objects, making it difficult to measure their effective 
temperatures T2 directly from observations, especially in the brown dwarfs r egime: in our sample , 
only two systems (HT Cas and DV UMa) have measured T2 available in the lRitter Kolbl ([2003 ) 
catalogue. One way to estimate T2 is to use the strong dependence of M2 on T2. Assuming that 
— M2 = (Mi) during the quiescence phase of DNe, we combine equations (|19j) and (123p to derive 

Ti pi/a \L 2 ) VO-9M y v ' 

Since donors with mass below the MHBM have L g /L 2 = 1, T 2 jT\ depends only on the orbital 
period and the effective mass-radius exponent. For masses above the MHBM, on the other hand, 
an accurate estimation of T 2 requires to consider the effect of (T g / L 2 )~ 1 ^ term, which can not be 
determined observationally. 



Using the method described in £|2.4[ we calculate L g /L 2 for the sample CVs listed in Table 
[2j For this, we take £ E = 0.21 ± 0.1 for systems with M 2 < 0.063 M and £ E = 0.64 ± 0.02 
for M2 > 0.063 M in accordance with the adopted empirical MRR. The bloating factor of each 
system is calculated relative to the 1-Gyr isochrone. We then calculate T 2 from equation (j24p and 
plot the results in Figure 0> as filled circles with errorbars; these are also tabulated in Column 
(9) in Table [2j It is remarkable that the estimated effective temperatures, except for SDSS 1507 
and DV UMa, are very close to the semi-empirical P-T 2 relation, suggesting that equation (I24p 
is in fact a good way to estimate T 2 when it cannot be measured observationally. A small orbital 
period of S DSS 1507 implies th at this system is most likely formed directly from a WD/brown 
dwarf pair ([Littlefair et al.l 120070 or a member of the old halo ([Patterson et al.ll2008l ). rather than 
a conventional DN. In the case of DV UMa, T 2 w 2400 K calculated from equation fl2U) is ~ 30% 
smaller than the observed e ffective temperature T 2 ~ 3200 ± 100 K, corresponding to a spectral 
type of M4.5±0.5, from the iRitter fc Kolbl (|2003l ) catalogue. Such a big difference in T 2 cannot 
be due to ambiguity in (£k + 1/3) 1 / 4 or (L g / 'L 2 )~ l / A terms since both are close to unity for this 
object. ITownslev fc Bildstenl (|2003l ) noted that the best-constrained quantity from a measured Ti 
is the accretion rate per unit area, (mi) = (Mi)/47ri?i, while the subsequent estimation of (Mi) 
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Fig. 7. — Comparison between our semi-analytic results (lines) and the 10 observed eclipsing CVs 
(filled circles) listed in Table [2] for (a) the mass-loss rate and (b) the effective temperature. The 
dashed and solid lines are for the Bl and B5 sequences, respectively, with the shaded and hatched 
regions denoting the uncertaities. In (6), Ti is estimated from equation (124p assuming that mass 
is conserved (i.e., — M2 = (Mi)). 
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from (riii) requires a certain MRR of WDs, which is poorly known for CVs. The WD mass in DV 
UMa is highest among the sa mple, and it is highly probable that the MRR for massive WDs differs 
fromiJi oc Mf 18 adopted bv Ebwnslev & Gansickel (j2009h . 



Based on the donor masses. iLittlefair et al.1 (|2008l ) claimed that SDSS 1035, SDSS 1501, SDSS 
1507, and SDSS 1433 possess post-period minimum CVs. Although period bouncers should have 
low M2 , the prec i se m ass of a particular secondary does not signify its evolutionary history or 
status ([Patterson! 12009 ) . The evolutionary status of such objects can be clarified more clearly 
by their positions, for instance, on the P-T2 or P-M2 diagram where the boome rang effect can 



be see n. Figure [7] suggests that SDSS 1501 and SDSS 1433 among the sample of ILittlefair et al 
(|2008l ) have not reached P m i n yet. This result does not depend on the stellar age, since these 
have — M2 > 10 -10 ' 5 M yr _1 for which the semi-empirical mass-loss rate is insensitive to the 
adopted isochrone. While the evolutionary status of SDSS 1035 is unclear, the fact that the P- 
M 2 dependence for the Bl sequence reconciles well with that due to gravitational wave emission 
indicates that it is likely a post-bounce CV. 



Kolb &: Baraffd (| 19991 ) argued that the tidal and the rotational distortions do not change the 



value of -Ptum appreciably because the internal structure of CV secondaries is only weakly affected 
by such perturbations. To show this, they compared the effective temperatures of the distorted 
and undistorted stars and found almost no change in T2. This result is totally expected since 
short-period CVs are fully convective stars. As we showed in §2.3[ the Hayashi theory gives that 
the reaction of T2 on the bloating of a donor is small for isolated, fully convective stars. We 
also show in Appendix |B] that the influence of an isentropic expansion on T2 caused by tidal and 
rotational distortions can be neglected in most practical cases. All of these suggest that the effective 
temperature does not reliably measure the amount of the stellar bloating, and consequently can 
not be used as an indicator of the degree of tidal and rotational distortions. 



5. Summary 

The mass-loss rate of donors is one of the most important parameters that govern the evolu- 
tion of short-period CVs. Since the mass-loss rate is not directly observable, most previous work 
employed either purely numerical methods or semi-empirical approaches that use the observed 
properties of CV accretors rather than donors. Yet, the confirmation of the resulting mass- loss 
rat es has to be m ade. In this work, we use the empirical MRR of short-period CVs constructed 



by iKniggd (|2006l ) as an input parameter and derive the semi-empirical relationship between the 
orbital period and mass-loss rate. We assume that the bloating of CV donors is caused mainly by 
the mass loss and associated thermal relaxation processes. We also calculate the responses of the 
effective temperature and luminosities on the bloating of the CV secondaries. We discuss our semi- 
empirical mass-loss rate in comparison with those from the previous numerical and semi-empirical 
studies, and in terms of estimating the effective temperature of the donors. Our main results are 
summarized as follows. 
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1. The semi-empirical P-M2 dependence shows that the mass-loss rate is essentially non- 
stationary. For P > 90 min, M2 is in the range of 10 -9,5 -10 -10 M yr _1 and decreases with 
decreasing P relatively slowly as M2 ot -P 3 - Below 90 minutes, on the other hand, M2 sharply 
declines towards P m m as a result of a strong dependence of M2 upon T%, emulating the P-T2 
relation. The P-dependence and amplitudes of the se mi-empirical M? are s ignificantly different 
from those predicted by the purely numerical models of iKolb &: Baraffd (|1999l ). The main reason 
for the discrepancies is small donor radii in the latter. Donors in short-period CVs are out of 
thermal equilibrium with r < 1. Thermal relaxation processes change the slope of the P-M2 and 
P—T2 diagrams significantly above P m m, making the boomerang shape less pronounced. Due to the 
sharp decline of T2 and M2 below 90 min, the observational detection of the boomerang shape on 
the P-M2 and P-T2 is highly unlikely. 

2. The semi-empirical P-M2 relation on t he bas i s of t he 1-Gyr isochrone is consistent with 
that estimated from the mass-transfer model of iRitterl (|1988l ) with the angular momentum carried 
dominantly by gravitational waves. Our semi-empirical P-M2 relation predicts that after emerging 
from the period gap, donors lose mass from 0.2 M Q to 0.075 M in ~ 1 Gyr, and reach the minimum 
period in ~ 1-2 Gyrs, suggesting that the ages of CV secondaries with 0.075 M© < M < 0.2 M© 
are less than 2 Gyrs. 

3. For 10 selected eclipsing CVs (listed in Table [2]) for which the reliable observational data are 
available, the accretion rates inferred from the observed effective temperatures of WDs are in good 
agreement with our semi-analytic mass-loss rates, indicating that our semi-empirical P-M2 relation 
is realistic. This also suggests that the effective temperature of a CV secondary can be estimated 
from the effective temperature of its WD primary through equation (|24p . When analyzed on the 
P-M2 or P-T2 plane, two CVs (SPSS 1501 and SDSS 1433) that were previously considered as 
post-bouncers bv lLittlefair et al.l (|2008l ) are most likely to be systems before reaching the minimum 
period. 
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government (MEST), No. 2009-0063616. 



A. Isentropic Expansion/Contraction of Polytropic Donors 



Even without considering the effect of mass loss, ISirotkin fc Kiml (|2009l ) showed that the tidal 
and rotational distortions alone make donors in CVs larger by about ~ 4-8% for n = 1.5 polytropes 
and ~ 2-3% for n = 3.5 polytropes, depending on the mass ratio. In this Appendix, we evaluate 
the enlargement factor for polytropes with arbitrary < n < 5. 



Consider a polytropic star with index n, mass M2, and pressure constant K. Its radius is i?2,i 
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in isolation. The MRR of the unperturbed polytropes is given by 

K = N nfi GM^- l)ln R% n)/n , (Al) 
where N n $ is the dimensionless pressure constant that varies weakly with n ( Chandrasekharlll933a ) . 



When it is situated in a CV, the presence of tidal and rotational perturbations causes it to expand 
or contract, depending on n, by adjusting the internal structure to restore hydrostatic equilibrium. 
Let the new equilibrium be achieved at the volume radius R2 with the corresponding dimensionless 
pressure constant N n . Since the global internal adjustment occurs in the dynamical timescale, 
shorter than the thermal timescale and the mass-loss time scale, one may assume that the specific 
entropy of the gas inside the donor does not change much under the influence of the distortional 
perturbations. 

For isentropic responses (i.e., K is unchanged), the enlargement factor is given simply by 

° isen = = ■ (A2) 



We use a self-consistent field method outlined in lSirotkin fc Kiml (|2009l ) to construct detailed inter- 



nal structure of polytropes in hydrostatic equilibrium. We then calculate N n fi and N n for polytropes 
in critical (i.e., Roche-lobe filling) configurations by varying n and the mass ratio q = M2/M1. Fig- 
ure [8] plots the resulting relative changes, Ai?2/-R2,o = «isen — 1, in the stellar size as a function 
of n for q = 0.1. The behavior of ai se n with M2 (or with the mass ratio q for M\ = 0.75 M ) for 
n = 3/2 polytropes is plotted in Figure [2] as a dotted line. Figure [8] shows that Ai?2 > for n < 3 
or n > 3.3, indicating that the perturbations make the donors bigger; the increase in the equatorial 
radius is larger than the decrease in the polar radius. On the other hand, critical polytropes with 
3 < n < 3.3 have AR2 < 0, since the perturbations decrease the polar radius much more than 
increasing the equatorial radius. 

Notice a discontinuity at n = 3 for which one cannot determine the stellar enlargement factor. 
For these polytropes, the central density is proportional to R% S and the ma ss becomes independen t 



of i?2) making it impossible to constrain the mass from its radius (e.g., IChandrasekharl Il933al ). 
Ai?2 diverges to either positive or negative infinity as n — > 3, indicating that a small change of n 
near 3 leads to a big difference in Ai?2- This suggest s that one shou l d be cautious when treating 



n = 3 polytropes numerically. Using SPH simulations. iRenvoize et al.l ([2002]) found that the critical 
configurations have an enlargement factor ~ 10-12% for n = 3 polytropes when q = 0.05-1.0. This 
is presumably due to some (unavoidable) numerical errors in the SPH simulations that somehow 
change their n = 3 polytropes ton~ 2.7-2.8 ones "effectively", resulting in Ai?2/^2,o = 0.10-0.12. 



B. Change of Effective Temperature due to Isentropic Expansion 



In Q2.3\ we describe how the effective temperature T2 of a fully-convective donor varies due to 
a homologous expansion using the Hayashi theory. Here we study how T2 responds to an isentropic 
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Fig. 8. — Changes in the volume radius of a polytrope due to the isentropic expansion under tidal 
and rotational perturbations against the polytropic index n, when the mass ratio is q = 0.1. A 
binary component expands when n < 3 or n > 3.3, while it contracts when 3 < n < 3.3. Note that 
n = 3 is a special case where the size of a polytrope cannot be constrained by its mass. 
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expansion. We assume that the pressure Pj in the whole stellar interiors below photosphere depends 
on temperature as 

_ ^/ 2T 5/ 2 

where [i is the mean molecular weight and 5i is the gas constant. The photosphere is determined 
at the radius where the optical depth is T p h = 2/3. If the radiative pressure is unimportant, the 
photospheric pressure, P p , is given by 

P P (r ph ) « (B2) 



where k is the opacity and g = GM2/R2 is the mean surface gravity. Using a simple opacity law 
of the form k = KoP5T b , we obtain 

^"(i^)* (B3> 

Since we assume that the donor is in hydrostatic equilibrium, the photospheric pressure should 
equal the interior pressure at the interface between the convective interiors and the radiative pho- 
tospheric layer. Assuming K = Kq, equations (IBlj) and (|B3j) are combined to yield 

T 2 



T21 



« ise r 26+5 . (B4) 



In contrast to the Hayashi theory, the effective temperature decreases due to the donor expansion. 
Nevertheless, T2 / /T2 i o = 0.99 for ai se n = 1.05, a = 0.5, and b = 4.5, typical values for CV secondaries, 
so that the change in T2 due to the tidal and rotational distortions are still very small and can thus 
be neglected in practice. 
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Table 1. The semi-empirical donor sequence for CVs 



l-Gyr 5-Gyr 



M 2 


R-2 


P 


t + t 


logZ/2 


T 2 


log 


Mi 


T 


t + t 


log Z-2 


T 2 


log 


M 2 


T 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


0.040 


0.099 


82.769 


3.50 


-4.625 


1273 


-11.10 


0.24 


13.08 


-5.513 


763 


-11.99 


0.24 


0.044 


0.101 


81.252 


3.10 


-4.435 


1407 


-10.94 


0.23 


9.74 


-5.317 


846 


-11.83 


0.23 


0.048 


0.103 


79.898 


2.80 


-4.293 


1513 


-10.83 


0.23 


7.31 


-5.167 


914 


-11.71 


0.23 


0.053 


0.105 


78.676 


2.55 


-4.194 


1588 


-10.76 


0.23 


5.47 


-5.011 


992 


-11.58 


0.23 


0.057 


0.107 


77.565 


2.33 


-4.121 


1645 


-10.73 


0.24 


4.08 


-4.880 


1062 


-11.49 


0.24 


0.063 


0.108 


76.021 


2.05 


-3.945 


1804 


-10.67 


0.29 


2.64 


-4.507 


1305 


-11.23 


0.29 


0.065 


0.111 


77.472 


1.95 


-3.828 


1905 


-10.61 


0.32 


2.30 


-4.316 


1439 


-11.09 


0.32 


0.069 


0.116 


79.730 


1.80 


-3.679 


2035 


-10.56 


0.40 


1.87 


-4.108 


1590 


-10.99 


0.40 


0.075 


0.121 


82.611 


1.70 


-3.481 


2226 


-10.39 


0.41 


1.72 


-3.649 


2021 


-10.55 


0.41 


0.078 


0.124 


84.036 


1.65 


-3.384 


2327 


-10.33 


0.44 


1.67 


-3.461 


2226 


-10.39 


0.43 


0.082 


0.128 


86.094 


1.59 


-3.255 


2466 


-10.21 


0.44 


1.60 


-3.282 


2427 


-10.23 


0.44 


0.086 


0.133 


88.096 


1.56 


-3.146 


2584 


-10.11 


0.45 


1.57 


-3.156 


2570 


-10.12 


0.44 


0.090 


0.137 


90.046 


1.53 


-3.063 


2671 


-10.04 


0.46 


1.54 


-3.064 


2668 


-10.04 


0.45 


0.095 


0.141 


91.948 


1.50 


-2.987 


2751 


-9. 


.98 


0.46 


1.52 


-2.986 


2751 


-9 


.98 


0.46 


0.099 


0.144 


93.805 


1.47 


-2.925 


2812 


-9. 


.91 


0.47 


1.49 


-2.924 


2813 


-9 


.93 


0.47 


0.103 


0.148 


95.619 


1.44 


-2.870 


2865 


-9. 


.90 


0.48 


1.46 


-2.869 


2866 


-9 


.89 


0.48 


0.107 


0.152 


97.394 


1.41 


-2.821 


2908 


-9. 


.86 


0.49 


1.43 


-2.820 


2910 


-9 


.86 


0.49 


0.111 


0.156 


99.132 


1.38 


-2.778 


2946 


-9. 


.84 


0.51 


1.40 


-2.777 


2948 


-9 


.83 


0.50 


0.116 


0.160 


100.835 


1.36 


-2.736 


2983 


-9. 


.81 


0.52 


1.38 


-2.735 


2984 


-9 


.81 


0.52 


0.120 


0.163 


102.504 


1.34 


-2.699 


3012 


-9. 


.79 


0.53 


1.36 


-2.698 


3014 


-9 


.79 


0.53 


0.124 


0.167 


104.142 


1.32 


-2.666 


3037 


-9 


.77 


0.54 


1.34 


-2.665 


3038 


-9 


.77 


0.54 


0.128 


0.170 


105.751 


1.30 


-2.636 


3057 


-9. 


.76 


0.56 


1.32 


-2.636 


3058 


-9 


.76 


0.56 


0.132 


0.174 


107.331 


1.28 


-2.605 


3080 


-9. 


.74 


0.57 


1.30 


-2.605 


3081 


-9 


.74 


0.57 


0.137 


0.177 


108.885 


1.26 


-2.575 


3104 


-9. 


.73 


0.58 


1.28 


-2.574 


3105 


-9 


.72 


0.58 


0.141 


0.181 


110.413 


1.24 


-2.547 


3124 


-9. 


.71 


0.59 


1.26 


-2.546 


3125 


-9 


.71 


0.59 


0.145 


0.184 


111.916 


1.22 


-2.521 


3141 


-9. 


.70 


0.61 


1.24 


-2.520 


3143 


-9 


.70 


0.60 


0.149 


0.188 


113.396 


1.20 


-2.497 


3155 


-9. 


.69 


0.62 


1.22 


-2.496 


3157 


-9 


.69 


0.62 


0.154 


0.191 


114.854 


1.18 


-2.472 


3174 


-9. 


.68 


0.64 


1.20 


-2.471 


3176 


-9 


.68 


0.63 


0.158 


0.194 


116.291 


1.16 


-2.447 


3191 


-9. 


.67 


0.65 


1.18 


-2.446 


3193 


-9 


.67 


0.65 


0.162 


0.198 


117.706 


1.15 


-2.425 


3206 


-9. 


.66 


0.67 


1.16 


-2.424 


3208 


-9 


.66 


0.66 


0.166 


0.201 


119.103 


1.13 


-2.403 


3219 


-9. 


.66 


0.68 


1.15 


-2.403 


3220 


-9 


.65 


0.68 


0.170 


0.204 


120.480 


1.11 


-2.383 


3230 


-9. 


.65 


0.70 


1.13 


-2.383 


3231 


-9 


.65 


0.70 


0.175 


0.207 


121.839 


1.10 


-2.365 


3240 


-9. 


.65 


0.72 


1.11 


-2.364 


3241 


-9 


.65 


0.72 


0.179 


0.211 


123.180 


1.08 


-2.345 


3252 


-9 


.64 


0.73 


1.09 


-2.344 


3253 


-9 


.61 


0.74 


0.183 


0.214 


124.505 


1.06 


-2.326 


3263 


-9. 


.63 


0.74 


1.07 


-2.325 


3264 


-9 


.64 


0.76 


0.187 


0.217 


125.813 


1.05 


-2.308 


3272 


-9 


.62 


0.75 


1.05 


-2.307 


3274 


-9 


.64 


0.78 


0.191 


0.220 


127.105 


1.03 


-2.291 


3281 


-9. 


.61 


0.75 


1.03 


-2.290 


3283 


-9 


.61 


0.81 


0.196 


0.223 


128.382 


1.01 


-2.275 


3289 


-9. 


.60 


0.76 


1.02 


-2.274 


3290 


-9 


.61 


0.84 


0.200 


0.226 


129.563 


1.00 


-2.260 


3295 


-9. 


.60 


0.78 


1.00 


-2.260 


3296 


-9 


.65 


0.88 



Note. — M 2 and i? 2 are the mass and radius of the donor in the solar units, respectively; P is the orbital period 
in minutes; t is the time in Gyr elapsed since the lower bound P g & p of the period gap, with to = 1 Gyr being the CV 
age at Pga P ; T 2 is the effective temperature of the donor in degrees Kelvin; L 2 is the luminosity of the donor in Lq; 
M 2 is the mass- loss rate in units of Mq yr" 1 ; r is the dimensionlcss parameter defined in equation {T}. 



Table 2. Data on Eclipsing Cataclysmic Variables (P < 130 min) with estimated effective temperature of WD 



System 
(1) 


P 

(2) 




Mi 
(3) 






9 

(4) 




Ti 

(5) 




M 2 
(6) 






ti 2 

(7) 




Ail 

(8) 




T 2 
(9) 


Ref 
(10) 


SDSS 1507 


DD.6119 





910 -0 


070 
07(1 





0625I 


0004 
0004 


11000I 1250 





nc7 +0 

057I 


004 
004 





0971°, 


003 
003 


— 10.78_ 


33 
34 


1 rTon + 360 

1730I 310 


3 


bDbb 1433 


78.1066 





8681° 


007 
007 





r\f c\r\ "t~U 

0690 _g 


0030 
0030 


12800I g50 





0601° 


003 
003 





mn+0 
109 -0 


003 
003 


— 10.431°, 


14 
15 


1 r>nn + 160 

1990Ii 60 


3 


bJJob loUl 


81.8513 





800l° 


030 
030 





(\n 7HTU 
06 ( 0_q 


0030 
0030 


12500_ 950 





053 n 


003 
003 





ma+0 
108_ 


004 

004 


— 10. 33_ 


19 
20 


20° -220 


3 


o T~\ oo 1 nor 

SDbb 1035 


82.0896 





940l° 


010 
010 





0550!° 


0020 
0020 


imnn+1000 

10100I 1000 





052_ 


002 
002 





108-0 


003 
003 


-10.981°, 


18 
20 


1370I 150 


3 


SDSS 1502 


84.8298 





8201° 


030 
030 





10901° 


0030 
0030 


123001°™ 





0901q 


004 
004 





1311° 


003 
003 


-10.401° 


19 
20 


2260l2?o 


3 


SDSS 0903 


85.0659 





960l° 


030 
030 





11701:" 


0030 
0030 


1300011°!° 





1121° 


004 
004 





1411° 


003 
003 


-10.581°, 


19 
20 


2470+ 740 
z4 ' u -430 


3 


SDSS 1227 


90.6610 





8101° 


030 
030 





11801°, 


0030 
0030 


150001^50 





0961°, 


004 
004 





1401° 


003 
003 


-10.031°, 


20 
21 


2690l|™ 


3 


OY Car 


90.8928 





8401° 


040 
040 





10201° 


0030 
0030 


150001^ 





086lg 


005 
005 





1331° 


003 
003 


-10.101° 


30 
33 


25501™° 


1,2,3 


HT Cas 


106.0560 





610l° 


040 
040 





15001° 


0100 
0100 


140001}°°° 





09ll° 


020 
020 





1501° 


013 
013 


-9.66I 24 


3000! 4 ™ 


1.2 


DV UMa 


128.2800 


1 


0411° 


024 
024 





15101° 


0010 
0010 


2000011™° 





1571° 


010 
010 





2041° 


016 
016 


a ' 3 ' -0.18 


23801 4 q° 


1.2 



Note. — P is the orbital period in minutes; Mi and M2 are the masses of the WD primary and the donor in the solar units, respectively; 
q = M2/M1 is the mass ratio; Ti and T2 are the effective temperatures of the WD and the donor in degr ees Kelvin, respectively ; R2 
is the radius of the d o nor in the solar units; M2 is the mass-loss rate in units of Mq yr — 1 ; References: (1) IPatterson et al.l j2005Tl ; (2) 
iTownslev fc Gansicke! fa00$\ ; (3) iLittlefair et al.l ||200S| 1. 



